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1 Introduction 

Decoherence [1-3] describes how a quantum system evolves into a state which most closely re- 
sembles a classical state through an interaction with the environment. As information about 
the system is lost to the environment as perceived by some observer, the entropy of the system 
increases. Recently quantitative studies have been done on decoherence for interacting quantum 
field theories, see Refs. [4-10]. These studies rely on the fact that neglecting observationally 
inaccessible, non- Gaussian correlators will give rise to an increase in Gaussian entropy of the 
state, which is then taken as a quantitative measure of decoherence. A useful application of this 
approach is to quantitatively study the decoherence of cosmological perturbations, i.e. cosmo- 
logical decoherence. During inflation quantum fluctuations of the inflaton field - and thus scalar 
cosmological perturbations - generically decohere because of the interactions generated by the 
inflationary potential as well as because of the interactions between cosmological perturbations 
and other fields present during inflation, and as a result become more and more classical [ll]. 1 
The ultimate goal of this project is to find out what the effect of cosmological decoherence on 
the cosmological microwave background (CMB) spectrum is for specific inflationary models, and 
thus to make a connection between inflationary models and late time cosmological observables. 

In this paper we focus on cosmological decoherence in the framework of stochastic infla- 
tion [13-15]. In stochastic inflation the system is represented by the super-Hubble (IR) modes 
of the inflaton field, as these modes are amplified during inflation and they are the ones that 
ultimately induce temperature fluctuations in the CMB. The sub-Hubble (UV) modes form the 
environment and act as a (white) noise for the long- wavelength modes. 

1 In Ref. [12] it is stated that even in the absence of interactions the quantum fluctuations become indistin- 
guishable from classical fluctuations during inflation. The authors call this "decoherence without decoherence". 
Technically, the state becomes extremely squeezed, but the phase space area remains constant (pure state). Our 
notion of decoherence differs, in that we associate decoherence with growth of the phase space area and conse- 
quently the entropy of the state. 
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Central in the stochastic inflation scenario is the Starobinsky equation 

4>(t, x) = ~V (<P(t, x)) + F(t, x). (1.1) 

The field eft is the coarse-grained part of the complete field, i.e. it contains only the long 
wavelength modes with k <C Ha. H = a/a (a is the scale factor) is the Hubble parameter in 
inflation and V(<f>) is a potential, e.g. V((f>) = A0 4 /4! for a self-interacting scalar field in a chaotic 
inflationary scenario. The noise term F originates from the short wavelength modes which keep 
crossing the Hubble radius due to the accelerated cosmic expansion. Interactions between UV 
modes are neglected, such that these modes are uncorrelated over time, thus the noise is white, 

(F(t,x)F(t',x)) = ^5(t-t'), (1.2) 

where (•) denotes a (statistical) average. Although the coarse grained field is a quantum field 
with creation and annihilation operators, in inflation it nevertheless commutes with its first 
time derivative (oc the canonical momentum). In this sense the field <f> may be called a classical 
stochastic field [15]. The Starobinsky equation is therefore a classical Langevin equation for the 
stochastic field (f>, in contrast to the quantum field equation for the complete quantum scalar 
field. From Eq. (1.1) a Fokker-Planck equation for the probability distribution of <j) can be 
derived, which is then used to calculate correlators for <j). 

Perhaps the greatest success of stochastic inflation is that it recovers the leading infrared 
logarithms {i.e. ln n (a) terms) of the field correlators at each order in perturbation theory [13, 15— 
19]. Originally it was believed that the time dependent separation between short and long 
wavelength modes would also lead to cosmological decoherence, even in a free field theory (see 
e.g. Refs. [20, 21]). However, Habib [22] showed that the free stochastic equations do not lead 
to any decoherence. 

Although the Starobinsky equation (1.1) is excellent for describing the leading order be- 
havior of the coincident field correlators, it fails when one is interested in the phase space of a 
state. Technically the Starobinsky equation cannot be derived from an action, thus no canonical 
momentum and Hamiltonian can be constructed. Moreover from the Starobinsky equation (1.1) 
no off-coincident two-point functions can be derived. This is especially troublesome for our 
approach to decoherence from the viewpoint of increased entropy due to incomplete knowledge 
of observationally inaccessible higher correlators [4-9] . The entropy of the system is in this case 
described by the Gaussian von Neumann entropy [4, 5, 8] 



A fc + 1 ^ M fc + 1 \ _ A fc - 1 ^ fA k - 1 



d 3 k 
(2^)3 

where A& is the phase space area (h = 1) for a Gaussian state centered at the origin, 



Ai 



<|0(*, £)| 2 ><M*, £)| 2 > " (l(4>(t, k)w;(t, k) + 7T*(t, k)<f>*(t, k))f 



(1.3) 



(1.4) 



where (f)(t, k) is a Fourier mode of the field 4>(t, x) that describes the system, and n(t, k) is its 
conjugate momentum. The phase space area is equal to 1 for a pure state, which consequently 
has zero entropy. For a mixed state the phase space area grows and entropy increases. Mixing 
occurs in general for an interacting quantum field theory, such as a self-interacting quartic 
potential for the inflaton field, used by Starobinsky and Yokoyama [14]. However, because the 
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canonical momentum is absent in the Starobinsky equation (1.1), the system will be in a 'pencil 
state' with zero phase space area. The entropy for such a state is not defined. Obviously we 
have to improve the Starobinsky equation in order to describe the growth of the phase space 
area and entropy. 

The solution is a full Hamiltonian (phase space) formulation of stochastic inflation, which 
was first (for a different purpose) used by Habib [22] for the case of free quantum fields. The 
idea is to introduce coarse-graining for both the scalar field and its conjugate momentum. As a 
result there will be two noise terms in the Hamilton equations for the IR field and its momentum 
(see also Ref. [23] for a more recent discussion). This is very much like a quantum particle 
which experiences non-commuting noise in both the position and momentum direction. For 
illustrative purposes we treat the simpler model of this quantum stochastic particle in section 2. 
We calculate the phase space area of a free particle moving under the influence of two types 
of noise: classical white noise and instant quantum noise. The latter is closely related to the 
Hamiltonian approach to stochastic inflation which we discuss in section 3 for a free scalar 
field. By calculating the phase space area we confirm that a pure state remains pure after 
the Starobinsky coarse graining. In section 4 we include a self-interacting potential for the 
coarse grained scalar field. We compare the stochastic Hamilton equations for a self-interacting 
field theory to the Starobinsky equation. A numerical method is proposed for solving the 
stochastic Hamilton equations which may present a first quantitative calculation of cosmological 
decoherence. Finally in appendices A and B we give basics of the quantisation procedure and 
the Wigner function approach to stochastic inflation, respectively. 



2 Classical white noise versus instant quantum noise 

In this section we study how the phase space area of a quantum particle changes under the 
influence of environmental noise by considering a simple toy model of a free quantum mechanical 
particle where both the position and momentum of the particle experience some environmental 
noise. These two noises induce kicks in the position and momentum direction. The Hamilton 
equations for a quantum particle moving under the influence of the noise terms are 

x = ^ + F 1 (t) (2.1) 
m 

p = -muj 2 x + F 2 {t). (2.2) 

The time dependent noise terms F\(t) and are environmental (white) noise, characterized 

by the noise correlators 

(Fi(t)) = (2.3) 
(F i (t)F j (t'))=f ij (t)5(t-t'), {i,j} = l,2. (2.4) 

The noise F2 is the ordinary noise term for a classical Brownian particle. Physically it can 
be viewed as a (big) particle that exchanges small bits of momentum by moving through an 
environment of smaller particles. The other noise term F\ has no classical analogue, as it 
represents a small infinitesimal displacement of the particle. However, it allows us to more 
easily draw an analogy with the stochastic scalar field in the Starobinsky picture. We will come 
back to this in the next section. By substituting Eq. (2.1) into Eq. (2.2) we find an equation 
for x 

x + uj 2 x = Fi + — F 2 . (2.5) 
m 
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We can now solve for x and p by considering the righthand side of (2.5) as a source term. The 
full solutions of Eqs. (2.1-2.2) are then 



x(t) = x (t) + / dt'G mt (t,t') 



d v F x {t') + -F 2 (t') 



m 



p{t)=p {t)+m / dt'd t G iet (t,t' 



d t ,F l {t') + -F 2 {t') 
m 



mFi(t), 



(2.6) 
(2.7) 



where xo(t) is the homogeneous solution of Eq. (2.5), po(t) 
function is 



mxo(i) and the retarded Green's 



G ret (t,t') = 9(t-t 



, sin[w(t - t')\ 



(2. 



So far we have not yet discussed canonical quantisation. The position x and momentum p are 
quantum operators which satisfy the canonical commutator [x,p] = i (setting H = 1). There are 
two ways to implement the quantum character of x and p into the solutions of the Hamilton 
equations (2.6-2.7): in Case 1 the homogeneous solutions xo and po contain the quantum nature 
of x and p. The noise terms F\ and F 2 are consequently classical white noise terms. In Case 2 
the quantumness of x and p originates from the noise terms, in the sense that F\ and F 2 do not 
commute, i.e [F\, F 2 ] ^ 0. In order to satisfy the canonical commutator the noise terms have to 
be of the special type of instant quantum noise. 

In both Case 1 and Case 2 we are interested in the phase space area A of the quantum 
particle, 

(x(t) 2 )(p(t) 2 ) - <^{*(t),p(t)}> 2 ] . (2.9) 

The phase space area can be used to calculate the Gaussian von Neumann entropy of the system 
(see for example Refs. [8, 24]), 



Srr = — — — In 



A + l 



A- 1 



In 



A-l 



(2.10) 



For a pure state the phase space area is minimal, i.e. A = 1, and consequently the entropy is 
zero. In the case of a free quantum mechanical particle without environmental noise the phase 
space area is minimal as there are no (environmental or internal) interactions that could lead 
to a growth of A. On the other hand, the phase space area in general increases once we include 
environmental noise. Our goal is to find out how the phase space area grows in both Case 1 
with classical white noise and Case 2 with instant quantum noise. We will now consider these 
two cases in more detail. 



2.1 Case 1: Classical white noise 

In the first case, the homogeneous solutions in Eqs. (2.6-2.7) are the usual quantum harmonic 
oscillator solutions expressed in terms of creation and annihilation operators., 



x (t) 
Po(t) 



1 



2mw 



moj 



(2.11) 
(2.12) 
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with [a, at] = 1 and [a, a] = [a\a^] = 0. These solutions satisfy the commutation relation 
[xo,po] = Provided that the noise terms F\ and F2 are classical, i.e. commuting, the full solu- 
tions x and p in (2.6-2.7) will satisfy the canonical commutator [x,p] = i. Now we can calculate 
the field correlators. The environmental noise and the free particle solutions are uncorrelated, 
i.e. (x Fi) = (p Fi) = 0. Using Eqs. (2.6-2.8) and Eq. (2.4) we find 

(x(tf) 
{Pit?) 



1 



2 TOW 

tow 



+ 



d£ 



ci/n(t') + s 2/22(0 



m 2 w 2 



+ sece 



(/12(f) + /21(f)) 



■ 2 2 

+ TO W 



dt' 



c 2f (J\ 1 -2 /22(f) 



{x(t),p(t}} y = mw / dt' 



, /,/N . /22(f) 

-cosefil{t ) + C6»S6> — + 



to^w 



to 2 w 2 ' 2 



(eg 



TOW 

(/12(f) + /21(f)) 
TOW 

2 , (/12(f) + /21(f)) 

9 e) 

TOW 



(2.13) 
(2.14) 

(2.15) 



where s# = sin0, eg = cos^ and 9 = ui(t — t'). We take the noise to be classical white noise. 
This means that the noise correlators in Eq. (2.4) are constant and that [Fi,^] = 0, i.e. 



/«(*) = /, 



13 ) 



fl2 - / 



21 



0. 



(classical white noise). 



(2.16) 



If we substitute these noise correlators (2.16) in the correlators Eqs. (2.13-2.15) we find the 
familiar result for ordinary Brownian motion that at late times (x(t) 2 ) oc t. Furthermore the 
variance of the momentum (p(t) 2 ) oc t. The phase space area (2.9) of the state becomes 



A = 1 + 2tow /n + 



22 



TO 2 W 2 



1 2 2 
+ TO W 



/ll + 



/22 



TO 2 W 2 



sin 2 cot 



to* 



122 



'11 



TO 2 W 2 



2 { (/12 + /21) 2 
TO 2 W 2 



At late times this implies 



A ~ 1 + TOW ( fn + 



/ 



22 



TO 2 W 2 



(2.17) 



(2.18) 



We clearly see that the phase space area of the initially pure state eventually grows linearly in 
time. As a consequence the late time entropy increases as S ~ 1n(t). These results are all well 
known for classical Brownian motion. The difference for the quantum Brownian particle is that 
the momentum as well as the position of the particle experiences a random environmental noise. 
As a result, the increase in phase space area depends equally on both noise correlators. 
To draw a connection with classical Brownian motion, we present here the original Langevin 
equation for classical Brownian motion, 

mx + (x = F(t), (2.19) 

where m is the mass of the particle and £ is a friction coefficient. The noise terms satisfy 

(F(t)} = 0, (F(t)F(t')} = 2(k B T5(t - f ), (2.20) 

with ks the Boltzmann constant and T the temperature of the particle's environment. The 
prefactor on the right-hand side of Eq. (2.20) is fixed by the fluctuation-dissipation relation. 
The classical Langevin equation (2.19) gives the familiar result that the displacement at late 
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times is proportional to y/i, specifically ((x — xq) 2 ) ~ 2Dt, where D = UbT/Q is the Einstein 
diffusion constant. On the other hand the variance of the momentum is (p 2 ) = m 2 (x 2 ) ~ m/c^T. 
This is to be contrasted with the variance of the momentum for our model (2.1-2.2), where we 
find from Eq. (2.14) that (p 2 ) oc t. The difference originates from the friction term in the 
classical Langevin equation (2.19) which limits the growth of (p 2 ). Technically the friction term 
is always present in a thermal bath because a fluctuation-dissipation relation must be satisfied 
for a Brownian particle. 

2.2 Case 2: Instant quantum noise 

In the second case the quantumness of x and p in the general solutions (2.6-2.7) originates from 
the noise terms F\ and F%. The homogeneous solutions are the usual classical solutions for a 
free particle 

xo(t) = x m cos[uj(t - t- m )] H — siu[o;(£ - t- m )] (2.21) 

mu 

Pin 

Po(t) = —mujx m sin[w(t — t- m )\ H cos[w(t — £i n )], (2.22) 

UJ 

where x; n and Vi n are the position and velocity of the particle at the initial time t{ n . Now, using 
the complete solutions (2.6-2.7) we find that the canonical commutator is 

[x(t),p(t)] = J'dt' [/ u (0 - / 21 (0] • (2.23) 

This suggest that our noise correlators (2.4) act only at one instant in time, and that the noise 
terms commute in a specific way, i.e. 

fij(t) = fijS(t — tj\r)j /l2 — /21 = h (instant quantum noise), (2.24) 

where t^ is the time at which the noise acts. Physically, one can view this is a classical particle 
moving freely, then being kicked by an infinitesimal sheet of smaller particles at time tjy. For 
classical Brownian motion with instant noise, the observer will see that the particle receives a 
small kick of momentum. On the other hand in the quantum mechanical case with noise in 
position and momentum (2.1-2.2), the particle experiences a kick of both its momentum and its 
position, i.e. a quantum kick. As a consequence of this kick, the phase space area of the particle 
becomes 2 

fu + /21 



A 2 = 46(t-t N ) 



fuf: 



22 



(2.25) 



Physically the result can be explained in the following way (see Fig. 1): at a time t < tjy 
the particle follows a classical trajectory in phase space described by the classical homogeneous 
solutions (2.21-2.22). The phase space area, expressed in units of h, is a typical quantum 
mechanical object. For a classical state the phase space area is zero, which can easily be seen 
by substituting the classical solutions (2.21-2.22) into Eq. (2.9). 



2 When (a;) 7^ and/or (p) 7^ 0, the correct definition of the phase space area is the following generalisation of 
(2.9): 

(Ax 2 )(Ap 2 }-(^{Ax,Ap}y 

where Ax = x — (x) and Ap = p — (p). Thus we must subtract any expectation value of x and p, which in this 
case are the classical homogeneous solutions xq and po in (2.21-2.22). 
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Figure 1: Trajectory of a particle influenced by instant quantum noise. From time ti n the particle 
follows the classical trajectory of a free particle in (2.21-2.22). At time ijv the quantum noise kicks 
in and the particle becomes aware of its quantumness. After the phase space area A is nonzero as 
described in (2.25). The classical trajectory is still shown as a dashed line. 



Now, at the time tjy the noise instantly kicks in. Through the noise the particle sud- 
denly becomes aware of its quantumness, and its phase space area becomes nonzero. In this 
sense the instant quantum noise provides a way of writing quantum mechanics as a stochas- 
tic theory. For a consistent quantum mechanical nature of the particle, we find the condition 



/ii/r- - ' 



122 - 2 



> 1, such that entropy is well defined. 



After the time t]y the phase space area will generally be nonzero. Consequently, the entropy 
increases by some amount due to the quantum kick of the particle. The Case 1 of quantum 
white noise means that the particle continuously receives quantum kicks. As we have seen, the 
phase space area grows linearly in time in that case. 



A special situation occurs when 4 



/11/22 - ( 



1. In that case the phase space 



area is minimal and the initially classical particle is now in a pure quantum mechanical state. 
As we will see in the next section, in stochastic inflation [13] the modes for a noninteracting field 
exhibit a similar behavior. In the phase space approach to stochastic inflation instant quantum 
noise is present in the Hamilton equations for the super- Hubble, or infrared (IR) modes of the 
inflaton field. The noise appears because of a special type of coarse graining, which is based 
on integrating out the sub-Hubble, or ultraviolet (UV) modes of the field. The properties of 
the noise are however extremely special, in the sense that the ultraviolet modes that cross the 
horizon are in a pure state, thus there is no entropy generation. We will examine this now in 
greater detail. 



3 Correlators and entropy for free coarse-grained scalar field 

We start with a field theory of a free massless scalar field in an expanding universe. Details of the 
theory and the quantized solutions of the fields and momenta can be found in Appendix A. In 
the phase space approach to stochastic inflation [22] the scalar field (A. 8) and momentum (A. 9) 
are separated into super- and sub-Hubble modes, or equivalently IR and UV modes, 

<& = 4> + ip (3.1) 
n = 7r + 7v (3.2) 
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The super-Hubble modes <j> and ir^ of the scalar field and momentum are those modes with 
k >C Ha. We separate these modes from the rest by introducing a time dependent cut-off 
9(eaH — k), where e < 1, and 

0(t, x) = f -0^6(eaH - k) {u(t, k)e^a{k) + u*(t, jk) e -^ 5f at(fe)} (3.3) 

7T^(t, x) = a 3 J -0^9{eaH - k) |n(t, k)j ls a(k) + u*(t, jfeje^a^jfe)} . (3.4) 

Similarly, the short wavelength expressions </3 and -zr^ are the other modes with k > eHa. The 
time dependent momentum cut-off in the short wavelength terms leads to a stochastic noise term 
in the long wavelength equations. We can easily see this by substituting Eqs. (3.1-3.2) into the 
Hamilton equations (A.4-A.5) and using the field equation for the mode functions (A. 11). We 
find 

4>= -V + Fi(t,x) (3.5) 
or 

^ = a 3 ^ + F 2 (t,x), (3.6) 
where the stochastic noise terms are 

F 1 (t,x) = J ^03 5(A:-ea^)(eaF)-{u(t,fc)e^a(fc) + ^(t,A:)e- i ^a t (A?)} (3.7) 

F 2 (t, x) = a 3 J S^K* - eaH){eaHy {«(t, k)j t3 a(%) + u*(t, k)e~ il3 a< (k)} . (3.8) 

where u(t, k) are the scalar field mode functions defined in (A.8-A.9), and a(k) and a\k) are 
the annihilation and creation operators (A. 10). Note that the stochastic Hamilton equations in 
the field theoretical case (3.5-3.6) are (almost) the same as those in the quantum mechanical 
case (2.1-2.2). At first it was believed that the coarse-graining of the free field could lead 
to decoherence due to the time dependent split in long and short wavelength modes [20, 21]. 
Habib [22] however argued that the stochastic inflation type of coarse-graining for a free field 
does not lead to decoherence. Here we will put these arguments on a firmer footing by performing 
an explicit calculation of (potential) decoherence for coarse-grained free scalar fields. We will 
do this in a general accelerating FLRW universe without assuming a specific choice of the scale 
factor. 

Our first goal is to show that coarse-graining through a time dependent cutoff does not lead 
to a growth of the phase space area (1.4) and Gaussian entropy (1.3) for the super-Hubble modes 
of the scalar field. In order to do so we will be interested in the momentum space correlators 
for and ir^ as they determine the phase space area for every mode k. In principle we can 
proceed along the same lines as the quantum mechanical particle with quantum white noise in 
section 2. The (free) scalar field is namely nothing more than an infinite sum of simple harmonic 
oscillators. Analogously to (2.5), by combining (3.5) and (3.6) and transforming into a Fourier 
space, we can write the field equations for the modes (j)(t, k) and ^(i, k): 

k 2 - 1 r - - 1 

<l>(t,k)+3H<j>(t,k) + -^</>(t,k) = -^ (a 3 Fi(t,A:))- + F 2 {t,k) (3.9) 
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and solve for <j>(t, k) using the method of Green's functions like in Eqs. (2.6) and (2.7). However, 
there is an important difference with respect to the quantum mechanical case. We obtained the 
stochastic Hamilton equations (3.5-3.6) by coarse graining the homogeneous solutions (A.8-A.9) 
and substituting the coarse grained field and momentum (3.3-3.4) back into the original Hamil- 
ton equations (A.4-A.5). It would therefore make no sense to add the homogeneous solution 
of Eq. (3.9) to the coarse-grained solution. Physically speaking the homogeneous solution is 
absent, as the IR quantum field only exists once the UV field has crossed the Hubble radius, 
which is described through the quantum noise terms (3.7-3.8) alone. 

Mathematically speaking, we could solve Eqs. (3.5-3.6) without knowing the origin of the 
equations, but with knowledge of the noise terms (3.7-3.8). Because the quantumness of (f> and 
7r^ is contained in the noise terms F\ and F%, we could still add a classical homogeneous solution. 
However, as we have also seen in section 2.2, the phase space area is defined in terms of quantum 
fluctuations on top of the classical solution. Therefore any classical solution does not appear in 
the phase space area and can for our purposes be neglected. Thus, the full (quantum) solutions 
for 4>(t, k) and n^t, k) are just 



<f>(t,k) = J dt'G ret (t,t',k) \ {a' 3 F 1 (t',k)y + F 2 (t',k)\ (3.10) 
7ty(t,fc) = a(tf fdt' [d t G ret (t,t',k)} [(a ,3 Fi(t',fc))- +F 2 (t',k)] - a(t) 3 iq(t, fc) , (3.11) 



where a' = a(t'). The retarded Green's function is given in Eq. (A. 15) of Appendix A. Of course, 
if one substitutes the actual noise terms (3.7-3.8) into the above solutions, we will recover again 
the coarse grained solutions in Eqs. (3.3-3.4). At this moment however we will keep using the 
solutions for (j)(t, k) and 7r^(t, k) in terms of the noise terms. Here we wish to calculate the phase 
space area in terms of the noise correlators. This allows us to draw an analogue with the 
quantum mechanical case in the previous section. First let us calculate the correlators for the 
noise terms in momentum space, which we can write in the following way, 



with 



(Fi(t, k)F*(t', k)) = f iS 8(t - t H )5(t - t') , (3.12) 

hi = \u(t H ,k)\ 2 (3.13) 

f22 = a 6 (t H )\u(t H ,k)\ 2 (3.14) 

(/12 + /21) = \a 3 (t H ) (u{t H , k)u*(t H , fc) + u(t H , k)u*(t H , fc)) • (3.15) 



In these expressions tn is the time of the mode split, i.e. the instant of time when k = eHa. 
From these expressions one can also see that the noise terms F\ and F 2 do not commute, since 

[Fx(t, %), FZtf, k)] = (/i2 - f2i)5(t - t H )5(t - t') = iS(t - t H )5(t - , (3.16) 

such that, strictly speaking, F\ and F 2 are operator valued noises. We have used here the 
Wronskian condition (A. 12). The crucial observation of the noise correlators in Eq. (3.12) is 
that they are precisely of the instant quantum noise type that we discussed in section 2.2, 
see Eq. (2.24). With this knowledge we can now calculate the phase space area (1.4) for the 
solutions (3.10-3.11). After various partial integrations and use of the Wronskian condition, we 
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find 



Al = 4S(t-t H ) 



r , , fl2 + hi 

J 11 J 22 



(3.17) 



Notice that does not depend on /12 — /21 and hence the operator character of the noises is 
not relevant for the calculation of the phase space area and thus neither for the calculation of 
the Gaussian entropy of the state. Not surprisingly, the phase space area of the IR modes (3.17) 
is precisely equal to the phase space area of the free particle influenced by instant quantum 
noise (2.25). Thus also here, once the UV modes cross the Hubble radius (i.e. at time tn), the 
quantumness of the IR modes kicks in through a nonzero phase space area. 

The question remains: how big is the phase space area for the super Hubble modes of an 
interacting scalar field under the influence of noise coming from the sub Hubble modes? Let us 
first consider the free scalar field theory from this section. In this case we know exactly what 
are the noise correlators from Eq. (3.12). We can therefore simply calculate the phase space 
area (3.17) of the coarse grained scalar field by inserting the noise correlators (3.12-3.15). By 
using the Wronskian condition (A. 12) twice we obtain the result 

A 2 k = 6(t - t H ) = 6{eaH - k) . (3.18) 

Thus for the super-Hubble modes with k < eHa the phase space area is always exactly 1, i.e. 
the system is in a pure state. The Gaussian von Neumann entropy (1.3) is therefore 



1< "p[^'"(^)-¥ i "(¥j|-' <-) 

This reflects the general fact that introducing a time dependent mode separation of a free scalar 
field by means of the Heaviside step function does not generate any entropy. Even though 
the Starobinsky coarse graining introduces noise terms for both the field and its conjugate 
momentum, the properties of the noise are such that, once the modes cross the Hubble radius, 
the phase space area is minimal. The reason is that the noise appears here not because of some 
coupling to an external environment, but because sub-Hubble modes are continuously crossing 
the Hubble radius. 

The same results for the phase space area (3.18) and entropy (3.19) can be, of course, 
obtained by using the coarse-grained solutions (3.3-3.4). We emphasize that these solutions are 
the exact solutions of the Hamilton equations with noise (A.4-A.5). If we use these solutions 
from the start, we trivially find 

(\(f>(t, k)\ 2 ) = 9(eaH - k)\u(t, k)\ 2 (3.20) 
(|7ty(t, k)\ 2 ) = 9{eaH - k)a 6 (t)\u{t, k)\ 2 (3.21) 

{<f>(t, jb),7rj(t,jfe)} ) = 9{eaH -k)a 3 (t)(u(t,k)u*(t,k) + u{t,k)u*(t,k^J . (3.22) 

Substituting these propagators into the phase space area (1.4) gives as result again Eq. (3.18). 
The results here are thus general: once we have a free field and produce noise terms by introduc- 
ing a time dependent cut-off for the modes, the entropy of the coarse-grained field will not grow. 
Here we considered a massless field, but even if we add a (time dependent) mass term, in prin- 
ciple we can still find normalized mode functions u{t, k) and we can express the coarse-grained 
fields as in Eqs. (3.3-3.4). In general adding or changing a mass term in the Hamilton equations 
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only changes the squeezing of a state, not its phase space area. Besides adding a mass term, we 
could also add a quartic potential A0 4 /4! for the coarse-grained fields. In a mean field approach 
we can write the corresponding Hamilton equations for the super-Hubble field and momentum 

as 

0= -^ + Fi(t,z) (3.23) 

= « 3 (5 " m ' " ^(*)>) + x) • (3.24) 

The term m 2 - \{<p 2 (t))/2 is a time dependent mass term, and hence in a mean field approx- 
imation the field remains a free field and the entropy does not grow. We emphasize that, in 
the original Starobinsky's approach, both interactions between the sub-Hubble fields as well as 
interactions between the sub-Hubble and super-Hubble fields are neglected. One can show [15] 
that, for a scalar field with an arbitrary potential in the stochastic formulation of Starobinsky, 
one correctly captures all leading order logarithms [ln(a)] n of coincident correlators to all orders 
in perturbation theory. The stochastic approach has been extended to Yukawa theory [16] and 
to scalar electrodynamics [17-19], but no stochastic formulation of quantum gravity is as yet 
known. 

An important question is then what can source decoherence within stochastic inflation. The 
simple answer is: one source of decoherence are interactions between the infrared fields treated 
beyond the mean field approximation. This source is present in our Hamiltonian formulation of 
stochastic inflation and it is the basis for the next section. There are, of course, other sources 
of decoherence in the full quantum theory, the most notable one are interactions between the 
sub-Hubble and super-Hubble fields. Including these would require a significant modification of 
the framework of stochastic inflation, and will not be considered in this work. Some aspects of 
these interactions are studied in Ref. [25]. However, Matacz (incorrectly) neglects the leading 
noise terms F\ and F2 (3.7-3.8) that are already present in the free theory (3.5-3.6). Finally, 
we remark that, based on the results of the study conducted in Ref. [7], we do not expect 
interactions between the sub-Hubble modes to generate any significant entropy. 

4 Correlators and entropy for a self-interacting scalar field 

Our main interest is to study decoherence of scalar cosmological perturbations during inflation. 
Here we model this problem by using a toy model (which captures some, but not all important 
features of cosmological perturbations): a self- interacting scalar field with a quartic potential 
V((j)) = A0 4 /4! in an accelerating universe. Since we are interested in the phase-space area, 
a full Hamiltonian analysis is necessary. The resulting classical stochastic equations are (c/. 
Eq. (3.24)), 

<P= 7 ^+F 1 (t,x) (4.1) 

^ = aV 2 + F 2 (t,f)-o 3 ^ 3 . (4.2) 

We emphasize that here Fi and F2 are are not the quantum noises given in Eqs. (3.7-3.8), but 
they are classical noise terms, which can be obtained from (3.7-3.8) by replacing the creation 
and annihilation operators by classical random variables as follows, 



a(k) — > a c (k) , a'(k) — > a*(k) , 
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where ct*(k) is simply the complex conjugate of a c (k). If the initial state is taken to be pure and 
free, then the classical (commuting) random variables a c (k) = M[a c (k)] + iQ[a c (fc)] are drawn 
from a Gaussian phase space distribution, which can correspond to a pure state (with the surface 
area A = 1), which can be squeezed and/or displaced, or to a mixed state with A& > 1 (at least 
for some k). This prescription resembles in spirit the methods used for studying preheating after 
inflation [26]. Several comments are now in order. 

Replacing the noises F\ 2 m (4.1-4.2) by their classical counterparts makes sense as long 
as correlators that involve anticommutators are large when compared with those that involve 
commutators of fields (and their momenta), and this is justified for the super-Hubble modes 
(see also Ref. [23] for a discussion about classicality) . The spirit of the classical approximation 
becomes even clearer when one uses perturbative two-particle irreducible (2PI) methods to 
study out-of equilibrium field dynamics. From e.g. Refs. [7] and [9] it follows that evolution 
of the statistical Green function (which involves a field anticommutator) depends also on (an 
integral over) the causal Green function (involving a field commutator). Now, if the causal Green 
function is much smaller than the statistical Green function, one can neglect the former in the 
equation for the statistical Green function, resulting in a classical approximation. This classical 
approximation is the 2PI analogue of the 1PI classical scheme advocated here. Ultimately, the 
validity of our scheme can be tested by comparing with the full (quantum) 2PI methods, which 
would require suitable adaptations of methods in Refs. [6], [7] and [9]. 

For a quantitative assessment of decoherence in cosmological settings one needs the Gaus- 
sian von Neumann entropy (1.3-1.4), which is a function of the phase space area A, and which in 
turn depends on correlators that involve only anticommutators. This then implies that solving 
the classical equations (4.1-4.2) should suffice. 

The main purpose of this work was to propose a simple method suitable for studying 
decoherence in cosmological settings. We postpone an actual numerical investigation to a future 
paper, where more emphasis will be given to the important question: how are cosmological 
observables influenced by the amount of decoherence of cosmological perturbations. 

Of course, the stochastic theory (4.1-4.2) must be discretized in order to solve these equa- 
tions numerically. The theory will be discretized in both position and momentum space at the 
level of the action in order to keep unitarity of the resulting (discrete) equations of motion (4.1- 
4.2). This is important, since such a discretization method guarantees that discretization will 
not induce a new source of decoherence. 

A further assumption in our stochastic approach is that there are no interactions between 
the UV and IR parts of the field. The validity of this approximation can be tested by making 
a comparison with a full 2PI evolution. This was firstly done by Giraud and Serreau [6] in 
Minkowsky space, and what remains to be done is to extend their work to inflationary spaces. 

Finally we note that alternatively one could calculate the Gaussian correlators in the phase 
space area from the Wigner distribution. This Wigner distribution can be solved from a Liouville 
equation which can be obtained from the stochastic Hamiltonian. For illustrative purposes we 
derive the Liouville equation for the toy model of the stochastic quantum particle in appendix B. 
Although an alternative to our proposed method, calculating cosmological decoherence through 
the Wigner distribution method is much harder due to the non-local nature of the self-interacting 
field in momentum space. 
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5 Conclusion and Discussion 



In this work we propose a generalised framework of stochastic inflation that is suitable for studies 
of decoherence of scalar cosmological perturbations during inflation. We take the Gaussian von 
Neumann entropy as a quantitative measure of decoherence. This means that the entropy is 
generated by neglecting observationally inaccessible higher order field correlators. 

In sections 2 and 3 free theories are discussed. Section 2 is mostly pedagogical. There we 
first consider the quantum stochastic particle with (white) noise in the direction of position and 
momentum. The canonical commutator can be satisfied by having a free quantum mechanical 
particle moving under the influence of classical white noise. In that case the phase space area 
increases linearly in time at late times. However, the quantumness of position and momentum 
can also originate from the noise terms themselves. The noise in that case is instant quantum 
noise, with non-commuting noise terms in the position and momentum direction which only act 
at one instant of time. In that case the initially classical particle receives a quantum kick and 
becomes aware of its quantumness by having a nonzero phase space area. The latter is closely 
related to the stochastic noise present in the phase space approach to stochastic inflation, of 
which we discuss the free case in section 3. The origin of the (quantum) noise is in fact that 
the sub-Hubble modes continuously cross the Hubble radius. It has the special property that it 
is unitary, i.e. once the modes exit the Hubble radius, the phase space area remains minimal 
and constant, implying that the modes are in a pure state. This special property guarantees 
that there is no decoherence for a free stochastic inflationary theory. In addition, there is no 
decoherence when interactions are treated in a mean field approximation. 

Next, in section 4, a classical Hamiltonian stochastic inflation for a self-interacting scalar 
field has been proposed, which should be suitable for studies of decoherence of cosmological 
perturbations during inflation. In our approach, the noises are taken to be classical stochastic 
commuting variables, such that the equations of motion are suitable for numerical treatment. 
The main advantage of our approach (when compared to more sophisticated 2PI methods) is 
its simplicity. Of course, it would be very nice to have an analytic argument concerning the 
accuracy of our method. In the absence of it, we point out that our approach to decoherence 
can be tested by the more sophisticated 2PI methods. 

Finally, we note that understanding cosmological decoherence is important, since a highly 
decohered state of scalar cosmological perturbations will shift the first acoustic peak and reduce 
the amplitude of the secondary peaks in the CMB radiation (see, for example, Ref. [27]). The 
ultimate goal of this project is to find out how strong the cosmological decoherence is in specific 
inflationary models, such as a chaotic A<3? 4 theory, and to find out what are the observable effects 
on the CMB. 
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A Free field quantization 

Here we recall the physics of canonical quantisation of a free scalar field in inflation. We start 
with the action of a free scalar field on a general background 



Sfree = J 'Xyf^g^ 'd^d^ . (A.l) 

We use the FLRW metric g^ u = diag(— 1, a 2 , a 2 , a 2 ) with a = a(t) the scale factor. The canonical 
momentum for is 

where C is the Lagrangian density, S = f d 4 xC. The Hamiltonian density becomes 

The free Hamiltonian H is related to the Hamiltonian density as H = f dPx'H. The Hamilton 
equations are 



* - as - > <*■"> 



Combining these two equations gives the free field equation 



<£> + - -^V 2 $ = 0. (A. 6) 

or 

We now quantize the system by imposing the canonical commutator 

[$(t,x),n$ (£,£')] = id 3 (x-x). (A.7) 

A general solution to the Hamilton equations (A. 4) and (A. 5) which satisfies this canonical 
commutation relation is 

m, x) = j ^3 ky l£ a(k) + u*(t, k)e- its a\k)^ (A.8) 

n*(t,£) = a 3 (t) J ^3 fc)e^a(fc) +u*(t, jfeje^a^fc)} . (A.9) 

The creation and annihilation operators satisfy the commutation relations 

[a{k),a ] {k')] = {2itf5 :i {k-k'), [a(k), a(k')] = , [a*(k), = . (A.10) 

The mode functions it(i, fc) obey the free field equation in momentum space 

k 2 

u + SHu + -rti = 0, (A.ll) 
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and satisfy the Wronskian normalization condition 



W[u(t, k),u*(t, k)] = u(t, k)u*(t, k) - u*(t, k)u(t, k) = (A.12) 

or 



The retarded Green's function for the field is 



G ret (i, x; t',x') = i9(t - 2), ^) 



i9(t - t') I -Z^j {«(*> k)u*(t', k) - u*(t, k)u(t', k)} . (A.13) 



which satisfies 



Q2 + 3F5t _ ly 2 ) G ret (t, x; t', f) = 6(t - t') 6 ^ X ' . (A.14) 

Its Fourier transform is 

G iet (t, t', k) = %6(t - t') {u(t, k)u*(t', k) - u*(t, k)u(t', k)} . (A.15) 

B Phase space area in stochastic quantum mechanics 

An alternative approach to the calculation of the phase space area and entropy of stochastic 
fields is to solve the quantum Liouville equation. Here we focus on a simple quantum mechanical 
example, but the approach can in principle be extended to coarse-grained inflaton fields. We 
will follow mostly the approach by Habib [22], which is in turn based on Kubo's approach [28]. 
We will clarify where necessary. 

We start with the stochastic quantum mechanical equations of motion, 

x=^ + F 1 (B.l) 

m 

p = -muj 2 x-V'{x) + F 2 , (B.2) 
which in principle follow from the Hamiltonian 

H = —p 2 + -muj 2 x 2 + V(x) + Fip - F 2 x. (B.3) 
2m 2 

x and p obey the canonical commutation relation [x,p] = ih. The time dependent noises satisfy 

(F l (t)) = (B.4) 
{F i (t 1 )F j (t 2 )) = f ij (t 1 )5(t 1 -t 2 ), {i,j} = 1,2. (B.5) 

The density operator p satisfies the von Neumann equation (the quantum Liouville equation) 

ihd tP =[H,p}. (B.6) 

In position space the density operator can be represented as 

/+oo r+oo 
dx \ dx'\x)p(x,x';t)(x'\. (B.7) 
-oo J — oo 
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The von Neumann equation projected onto position space becomes 

ihd t p(x,x';t) = |-^(^ 2 - d 2 ,) + ^(x 2 - x' 2 ) + V(x) - V(x') 

- iHF 1 (d x + 8 X ,) - F 2 {x - x')\p(x, x'; t) (B.8) 

where we have made used that p\x) = —ikd x \x) using the canonical commutator [x,p] = ih. 
Now we perform a Wigner transform of the density matrix. This gives the Wigner function 

W(Q,P;t) = f d^xe-^ h p[q+^, q -^-t^. (B.9) 
The new coordinates are defined as 

x + x' 

q -- 



2 

Ax = x — x . 

Substituting these coordinates into the von Neumann equation (B.8) and using that 

(Ax Ax \ 
q + ~y^~ ~Y' t ) 

Axe -ip*x/h = ih U-ip&x/h\ ^ 

we find that the equation for the Wigner function (B.9) becomes 

\d t + — d q - muj 2 qd p + ^-AV + 2F 1 d q + F 2 d p \ W(q,p;t) = 0, 
[ m in J 

where 

= V (q + jd^j -V (q-jdp) ■ 
We write Eq. (B.14) schematically as 

{d t + L + Lj + L N (t)} W(q,p; t) = 0, 



where 



Lq = —d a — mu) 2 qd p 
m 

Lj = ^AV 

in 

L N (t) = 2F 1 (t)8 q + F 2 (t)d 1 



V 



For specific potentials 



V=^^L I , l = ^q 2 d p -\h 2 dl > 

A 4 A 



V = '-x* — ► L I)2 = ^(h 2 qd 3 p - 4q 3 d p 
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We now strip off the time dependence in Eq. (B.16) by denning 

<r(q,p;t) = e {t - to){Lo+Ll) W(q,p;t), 

which leads to the equation 

Q t(J = _ e (t-to)(Lo+L I ) LN ^ e -(t-t )(L +L 1 ) a 
= Q(t)a 



(B.22) 



(B.23) 
(B.24) 



The formal solution of Eq. (B.16) is 



v(q,P;t) = T < exp 



to 



o-(q,p;t ) 



1+ f dt x n(ti)+ I dh f dt 2 n(t 1 )n(t 2 ) + ..^]a(q,p;t ), (B.25) 
Jt Jt Jto J 



where 



cr(q,p;t Q ) = W(q,p;t ), 



(B.26) 



and T denotes the time ordering operation. Now we want to average over the noise in Eq. 
(B.31) by using the noise correlators (B.5). We assume that the initial Wigner distribution is 
uncorrelated with the noise, 



(Fi(t)W(q,p;t )) = = {Fi{t)a{q,p;t Q )). 



We then find that 



(a(q,p;t)) = ( T <^ exp 



dtiQ(ti 



to 



{a(q,p; t )}. 



The exponential can be expanded into powers of O(t). We use that 

<«(*)> = o 

(n(ti)n(fe)) = s(h -t 2 ) e (*i-'o)( i o+L/) jD 2^ i)e -(t 1 -t )(L +L,)_ 



(B.27) 



(B.28) 



(B.29) 



where 

D 2 {ti) = 4/u(ti)^ + 2(/i 2 (ti) + f2i{ti))d q d p + f 22 (ti)d 2 p . (B.30) 

Without loss of generality we now set to = 0. Furthermore we use some shorthand notation 
where 



dtin(h), 



t r ti 



1"2 



dti / dt 2 0,(t 1 )n(t 2 ), etc. 
I! Jo 
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We therefore find that 



(a(q,p;t)) = ( T <^ exp 



+- 



t r t 



/ (Oi^ 2 )+/ / / / (niQ 2 )<Q 3 ^4> + .--J <<7(g,p;0)> 

= ^1 + 1 J' e {L 0+ L 1 )t 1D 2 {ti)e -{L 0+ L I )t 1 
l e (Lo+i/)ti D 2^ i ) e -(Lo+L / )t le (Lo+L / )t 2jD 2^ 2 ) e -(L +Li)t2 + ... | (a(q,p; 0)). 



In the third line we have performed Wick contractions of the four- and higher-point correla- 
tors. Notice that contractions other than the time ordered one vanish. The final line can be 
reexponentiated to the simple form 



(a(q,p;t)) = T ^ exp 



1 f dt 1 ^ Lo+L ^D 2 {t 1 )e^ Lo+Ll ^ 

2 Jo 



(a(q,p;0)). 



Differentiation of the last line gives 



d t (a(q,p;t)) = ^ +L ^D 2 (t)e~^ +L ^(a(q,p;t)). 



(B.31) 



(B.32) 



As Lq and Lj only depend on the classical quantities q and p, they can be taken out of the 
noise average of a. This allows us to reexpress Eq. (B.32) in terms of the Wigner function W 
by using Eq. (B.22), 



d t + L Q + L I -^D 2 {t)}\Y(< t . ! Kl)-i). 



(B.33) 



This is the generalised Fokker-Planck equation in phase space variables. To remind the reader, 
Lq, Li and D 2 are defined in Eqs. (B.17), (B.18) and (B.30), respectively. The Wigner distri- 
bution can be used to calculate field expectation values for operators of x and p, 



(0(x,p))= / dqdpO(q,p)W(q,p;t). 



(B.34) 



We emphasize here that the x and p on the lefthand side are operators, whereas q and p on the 
righthand side are phase space variables. One has to be careful when calculating phase space 
correlators for which the ordering of x and p is important. The prescription is to move all the 
p's to the right in the operator 0(x,p), and then apply Eq. (B.34). This is of no concern for our 
object of interest, the phase space area (2.9), which does not contain any correlators involving 
the canonical commutator. 

In principle one could make an ansatz for W(q,p;t) which solves the Fokker-Planck equation 
(B.33). In the case of a free particle the solutions for the Gaussian Wigner distribution and 
related correlators can be obtained analytically, see for example [8] and [24]. For an interacting 
theory the non-Gaussian Wigner distribution cannot be solved exactly and one has to resort to 
perturbative methods. 
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We note that the above Wigner formalism for a quantum mechanical particle can be ex- 
tended to quantum fields. This has been done by Habib [22] for free fields. In principle the 
Wigner distribution W and the phase space area can be calculated numerically in a discretized 
version of Eq. (B.33) for interacting quantum fields. However, we warn the reader that it would 
be very challenging to find the Wigner function even numerically due to the non-local nature 
of the interacting potential in Wigner space. Thus, although Eq. (B.33) provides an alternative 
method, our proposed method in section 4 of solving the stochastic Hamilton equations (4.1-4.2) 
numerically provides a simpler way to model cosmological decoherence. 
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